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ABSTRACT 

The quark-antiquark interaction from the NRQCD Lagrangian is studied in 
the Wilson loop formalism. 



1. The NRQCD Lagrangian 

Large radii quarkonia (i.e. excited heavy mesons like bb and cc) cannot be de- 
scribed in terms of perturbative QCD with the addition of leading nonperturbative 
effects encoded into local condensates.u This happens since in this situation the glu- 
onic correlation length cannot be considered infinitely large with respect to the other 
scales of the system. A successful description needs therefore a systematic inclusion of 
non-local condensates. A solution is provided by the so-called Wilson loop formalism. 
The non-local quantities are field strength insertions in the Wilson loop made up by 
the quark trajectories.&ta 

While a full relativistic QCD formulation in this formalism is still lacking^ such a 
formulation is possible for heavy quark bound states which are essentially described by 
non relativistic dynamics. These systems are characterized by a dynamical adimensional 
parameter, the quark velocity v , which is small and allows a classification of the energy 
scales of the problem in hard (~ m), soft (~ mv ) and ultrasoft (~ mv 2 ). Moreover, this 
provides a power counting scheme for the operators in the Lagrangian. The relation 
between v and the QCD parameters is unknown (for infinitely heavy quarks v coincides 
with a s , for realistic quarks it is the result of perturbative and nonperturbative effects) 
but irrelevant once the power counting scheme works. Due to the heavy quark mass m, 
at a scale fi between m and mv the physics is still dominated by perturbative effects. 
Therefore, in order to describe heavy quark bound states, it is possible to substitute 
the QCD Lagrangian with an effective non relativistic Lagrangian via perturbative 
matching at that scale. The new Lagrangian is simpler, since the hard degrees of 
freedom have been integrated out explicitly, but equivalent to the QCD one at a given 
order in a s and v. This effective Lagrangian is known as the NRQCD Lagrangian. a At 
order 0{y A ) the NRQCD Lagrangian ^describing a bound state between a quark of mass 
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where are the heavy quark fields. The coefficients evaluated at the 

matching scale \i for a particle of mass mj. They encode the ultraviolet regime of 
QCD order by order in a s . The explicit expressions and a numerical discussion can be 



found inB The power counting rules^l for the operators of Eq. ([I]) are Q ~ (mv ) 3 / 2 , 
D ~ mv, gA ~ mv 2 , gA ~ mv 3 , gE ~ m 2 v 3 and gB ~ w?v A . Four quark operators 
which are apparently of order t> 3 are actually suppressed by additional powers in a s 
in the matching coefficients and the octet contributions by an additional power in v 2 
on singlet states. Therefore in the following we will neglect these contributions with 
the exception of a term which mixes under RG transformation with the chromomag- 
netic operator contribution to the spin-spin potential.Ej We will call the corresponding 
matching coefficient d. 

2. The Wilson Loop Formalism 

The use of the Wilson loop formalism on the Lagrangian ([!]) was first suggested 
in.El Let us sketch the derivation of the heavy quark potential. The 4-point gauge 
invariant Green function G associated with the Lagrangian (P is defined as 

G(x l , y u x 2 , y 2 ) = (0\Q f 2 (x 2 )(j)(x 2 , x 1 )Q 1 (x 1 )Q\(y 1 )<f>(yi, y 2 )Q2(2/2)|0), 

where (p(x 2 ,Xi) = exp j— ig J ds (x 2 — x^A^xi + s(x 2 — £1)) j is a Schwinger line 

added to ensure gauge invariance. After integrating out the heavy quark fields, G can 
be expressed as a quantum-mechanical path integral over the quark trajectories^: 
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where the bracket means the Yang-Mills average over the gauge fields, F is the Wilson 
loop made up by the quark trajectories Z\ and z 2 and the endpoints Schwinger strings 



and 2/° = 2/i = — x\ = x\ = T jl. Reparameterization invariance0 fixes c 2 = C4 = 1. 
Assuming that the limit exists, we define the heavy quark-antiquark potential V as 

rT/2 

lim / dt U. For a discussion on the existence of a quark-antiquark potential we refer 

r -too | jLpT/2 

toHM Expanding in v (following the power counting given in the previous section) or 
in the inverse of the mass we get 
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VK(r) = P exp i^—ig(j> dz^A^z)^ is the non-static Wilson loop. The expansion of it 

around the static Wilson loop VT(r ) (r is a r x T rectangle) gives the static potential 
Vq = lim ilog(W(To)) /T plus velocity (non-spin) dependent terms.H and are 

T — >oo 

the spin and orbital angular momentum operators of the particle j. The matching 
coefficients are defined as 2c FS = c FS ± c F g. The "potentials" V±, V%, ... are scale 
dependent gauge field averages of electric and magnetic field strength insertions in the 
static Wilson loop: 
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where (( )) = ( W(T ))/(W(T )). 



3. Conclusions 



Two comments in order to conclude. The explicit expression for the potential 
(0) in terms of field strength insertions in a static Wilson loop is suitable for direct 
latticed and analyticE-3 evaluations. This makes the obtained result of particular impor- 
tance. Different vacuum models can be easily compared on the heavy quark potential 
predictions, once the Wilson loop average has been evaluated. 

The 0(v 2 ) leading order NRQCD Lagrangian L = Q\(iD + d 2 /2m 1 )Q 1 + anti- 
quark part does not contribute to (0) only with a static potential (with the exception of 
the perturbative contribution if evaluated in Coulomb gauge). Since the corresponding 

Wilson loop P exp j -ig £ dz°MA } » a function of the non-static loop r, its expan- 

sion produces in general velocity dependent terms as well.00'0 This is not surprising 
since the power counting scheme of NRQCD has to be considered as a leading order 
power counting scheme. An exact value in v cannot be assigned to each term of the 
effective Lagrangian at least before the soft and ultrasoft degrees of freedom have been 
disentangledJiHl 

An extensive analysis of the topic discussed here can be found in,ta with particular 
emphasis on the relevance of the matching procedure in order to have a consistent 
potential in the perturbative regime. 
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